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Abstract 

We study dynamical fermion effects in lattice QCD at finite tem- 
perature. The method adopted is basically the extrapolation from neg- 
ative flavour numbers already tested at zero temperature and based 
on the simulation of local bosonic theories, with an essential differ- 
ence. With an appropriate choice of the boundary conditions on the 
bosonic fields, called "bermions", it is possible to separate the 
breaking contribution of fermion loops to the effective action from the 
one conserving Z3: the former is simulated exactly at a fixed posi- 
tive and even flavour number, while the extrapolation from negative 
flavour numbers is made only on the Zj, invariant part of the action. 
We test this approach by comparing our results on a 16^ x 2 lattice 
with those from a hopping parameter expansion and our results on a 
16^ X 4 lattice with those of direct Monte Carlo simulations including 
the fermion determinant. 
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1 Introduction 



A new approach to the problem of estimating dynamical fermion effects in 
lattice QCD based on the extrapolation from theories with a negative num- 
ber of dynamical flavours was recently proposed and applied in numerical 
simulations of the theory at zero temperature. In that case the quenched 
approximation appears to reproduce most of the experimental pattern of the 
hadronic spectrum and unquenching corrections are small when the compar- 
ison with the quenched case is performed after a suitable shift of the bare 
parameters i. e. at fixed values of the lattice spacing and of the renormalized 
quark mass. The "bermion" method exploits the smoothness of dynamical 
flavour dependence studied at fixed renormalized quantities to make esti- 
mates in full QCD which are extrapolations of results obtained at negative 
flavour numbers, where fermions are replaced by "bermions", i.e. bosons 
with a fermion action [|l], ^ . 

In this paper we extend this approach to QCD at finite temperature. 
The unquenching effects are not expected in this case to reduce to a simple 
renormalization with some small residual corrections. The SU (3) pure gauge 
theory is invariant under the discrete group and the deconfining transition 
in the theory is associated with the spontaneous breaking of this symmetry. 
In full QCD the Z3 symmetry is broken explicitely by fermion loops and 
the corresponding phase transition could disappear for intermediate values 
of quark masses. The straightforward application of the bermion method 
to finite temperature QCD would require an extrapolation from negative to 
positive flavour numbers rif also for the breaking fermion loop contribu- 
tion to the effective action. This extrapolation seems quite unsafe because 
thermodynamical QCD observables, e. g. the critical temperature, are not 
expected to behave smoothly around Uf = 0. However it is possible, using 
appropriate boundary conditions on the bermion fields, to simulate theories 
where the Z^ breaking and the Z^ conserving contributions of fermion loops 
to the effective action get different weights corresponding to different "ef- 
fective flavour numbers". As a result, the Z3 breaking interaction can be 
simulated exactly for even and positive flavour numbers, while the Z3 con- 
serving dynamical fermion contribution to the action, which is expected to 
give mainly a renormalization effect, is extrapolated from negative flavour 
numbers. 

In section 2 we present the bermion action and the choice of the bound- 
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ary conditions which leads to the correct breaking interaction. In section 
3 we compare our results on a 16'^ x 2 lattice with those of a hopping param- 
eter expansion and our results on a 16^ x 4 lattice with those of direct 
Monte Carlo simulations including the fermion determinant for A^^ = 4 0]. 

2 The role of the boundary conditions 

The action for lattice QCD with Uf flavours of Wilson fermion is 

s[u,^,4^] = 5G[f/] + EE^.(^)75[Q^i](x) (1) 

j=l X 

where S'G[f/] is the gauge action, j is a flavour index and 'j^Q is the Dirac 
operator: 

1 1 3 

[Qipi]{x) = — 75^i(a;) - -75 U^,{x){l - -f^,)iJi{x + /i) 

^ ,1=0 

-^75 E Ulix - + - li) (2) 

After integration over the quark fields one obtains an effective action: 

S,jf[U] = Sg[U] - n;Tr(log75Q) (3) 

where the second term is a sum of one loop diagrams. These can be divided 
into two classes: to the first belong the loops which lie entirely inside the 
lattice or wrap around the edges of the lattice a multiple of three times 
closing through the boundary conditions and to the second those which wrap 
a number of times different from multiples of three. At zero temperature all 
winding loops are just a finite volume effect, but at finite temperature, i.e. 
at finite physical extent of the "time" direction, the loops winding in this 
direction and belonging to the second class do influence the existence of 
the deconfining phase transition and an extrapolation from negative flavour 
numbers cannot be expected to be smooth. However, there are choices of the 
boundary conditions which affect only the second class and which reproduce 
the correct fermionic Z3 breaking contribution also in the case where the 
fermion determinant is replaced by a bermion one. 
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The action for QCD with a negative flavour number is written in 
terms of nb = \nf\/2 bermions fields which are commuting spinors: 

rib 

= Sg[U] + E<^](^)Q(^'^)Q(^'2/)<^i(2/) (4) 
x,y,z j=i 

where j is a flavour index. We introduce a set of flavour dependent boundary 
conditions parametrized by the elements of the center of the gauge group Z3: 

t + Nt) = -Zj(j)j{x, t) Zj e Z3 (5) 

The contribution of rib flavours to the second class of bermion loops in the 
effective action can be expressed as a sum of terms which are characterized 
by winding number modulo three i of Polyakov loop-like terms Pi in the 
"time" direction. The full expression for Sejf can be parametrized as: 

"6 rif, 

S,ff = SG + n,Re[Po] + Re[{Y.zj)P^]+Re[{j2^])P2\ (6) 

where zj is the element of which fixes the boundary conditions of the j-th 
flavour. The property Zj=l implies that the maximum power of zj that 
can appear in eq. ^ is two. The expression above should be compared with 
the corresponding fermion case with n/ flavours and ordinary antiperiodic 
boundary conditions: 

Seff = Sg- Y[i?e(Po) + Re{Pi) + Re{P2)] (7) 

The bermion boundary conditions can change the signs of Pi and P2 
terms in the effective action: for example with two bermion flavours and the 
choice 

zi = exp(+27ri/3) 
Z2 = exp(— 27ri/3) 

one gets: 

Seff = Sg + 2i?e(Po) - i?e(Pi) - Pe(P2) (8) 

i.e. the contribution of two bermion flavours has been changed, for the second 
class of diagrams, into the one of two fermion flavours. We introduce a 
"normal" fermion number Uf and two "magnetic" ones, rrifi and mj2 to 
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distinguish between loops with winding number 1 and 2 respectively and we 
write the effective action in the form: 



Seff = Sa- f i?e(Po) - ^ReiPi) - -fReiP^) (9) 

In table 1 we give, for various numbers of bermion flavours, the values 
of Uf, rrifi and that are obtained with specific choices of the boundary 
conditions. The correct fermion magnetic flavour numbers can be obtained 
only in some special cases. The extrapolations to positive nj at fixed val- 
ues of rrifi and mj2 can be made only for a few values: for example, for 
rrifi = mf2 = 2 the available nf are —5 and —2. For the other cases the 
table includes a further choice of boundary conditions: the periodic ones 
{zj = —1). By adopting these conditions for some of the bermion flavours 
we can simulate theories where the values of m/i and m/2 are different and 
also the contributions of terms with winding number larger than 2 are mod- 
ified. In the case of moderately heavy quark masses where winding numbers 
higher than 2 can be neglected, the m/2 contributions can be seen as a small 
correction to the rufi terms. In this case thay can be extrapolated together 
with Uf terms. In the heavy quark mass case also the terms with winding 
numbers 2 can be neglected and the fate of m/2 terms ignored. 

In the coming section we will investigate the role of the winding number 
2 term: at Nt = 4 and for k > 0.14 we find some evidence for the relevance 
of such a term. 



3 Comparison with existing results 

We have tested our method at A^^ = 2 with the results of the hopping parame- 
ter expansion. According to ref. 0, the effective action in this approximation 
and at leading order reads: 

Nt 

S,ff = SG + Hj2ReTrY[ U^{x, t) + o{k^) (10) 

X t=l 

where H = 2mji(2/t)^*. The shift of the critical f3 value Pc{H) from (3c{0) 
can be predicted [Q in term of the jumps of Polyakov loop and plaquette 
expectation values in the pure gauge theory: 

(3,{H) = (3,{0) - (4.94 ± 0.75)H (11) 
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To leeding order in k only terms with winding number 1 contribute 
and only the value of mji is relevant. Periodic boundary conditions instead 
of antiperiodic ones are sufficient to promote bermion effects into fermion 
effects. In order to test the expression above, we have run different values of 
Hf (from —2 to —10), mji (from 2 to 6) and k (up to 0.071), while keeping 
the coefficient SrufiK^ of the Z3 breaking term in eq. |iy equal to a fixed 
value ranging from 0.01 to 0.08 and we have measured Polyakov loops and 
plaquettes. For this small values of k the effects of different values of n/, 
which are of higher order in the hopping expansion, should be negligible and 
in fact the values of the observables at fixed /? depend only upon Sm/iK^. 
We have determined (3c{H) by monitoring the jump in the Polyakov loop 
expectation value. The results, shown in fig. 1, are in remarkable agreement 



with the prediction of eq. |TT|, confirming the validity of the method in the 
hopping expansion regime. 

Increasing A^^ and going to higher values of k rapresent a more severe 
test for the method. We have compared our results for Nt = 4 with those of 
hybrid Monte Carlo simulations of ref. [^], and in particular we have studied 
the behaviour of the Polyakov loops and of the plaquettes around the phase 
transition point. In this case the rufi term starts at the order k^, the renor- 
malization effects due to "normal" fermion loops are not negligible anymore 
and we have to perform an extrapolation from negative values, at fixed val- 
ues of rrifi. As already discussed, the values of m/2 cannot be fixed at the 
fermionic ones at all values of n/. In figures 2-4 we present the behaviour as 
a function of j3 of the space-space plaquette 1-Pa = -^Re(Tr{UiU2Ulu\)) 
and of the absolute value of the Polyakov loop L, defined by: 

for different values of n/ and constant values ofm/i at k = 0.12, 0.14 and 0.16 
respectively, together with the extrapolations to full QCD with two fermions 
[uf = rufi = mf2 = 2). In figures 5 we present for the case k = 0.14 the 
results for the difference between space-space and space-time plaquettes. 

The procedure adopted for the extrapolation is the following: we keep 
fixed the values of k and of a given observable {P„, \P„ — Pt-\ or |L|) and 
extrapolate the corresponding /? as a function of n/. From previous applica- 
tions of the bermion method to QCD at zero temperature it is known that, in 
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the intermediate quark mass region, a good extrapolation in n/ is obtained 
by keeping fixed renormalized quantities. This is true in principle also for 
the simulations at finite temperature, however it must be noticed that our 
results refer, at least across the phase transition jump, to a rather heavy 
quark mass region where the renormalization effects are not so large to re- 
quire a non perturbative treatment and therefore the extrapolation at fixed 
bare parameters gives a good first guess. 

The results refer to the cases a, b, c and e of table 1 which correspond 
to m/i = 2 and four different negative values oinj. For small values of n the 
effect of the terms with winding number larger than one should be negligible 
and therefore, indipendently from the value of m/2 the data should align as 
a function of n/, as is the case for k = 0.12 and 0.14. Even if they do not 
align we can consider two possible extrapolations which take into account 
the relevance of terms. The first one is from the cases b and e where 
all the breaking terms are fixed at the fermion value rrifi = mf2 = 2. 
The other one is to extrapolate from cases a and c: it keeps fixed at the 
fermion value only the dominant Z3 breaking term {rrifi = 2), while mj2 is 
extrapolated from negative to positive values together with n/. The results 
of the two extrapolations are given in the figures: they agree with each other 
and with the results of direct hybrid Monte Carlo simulations Q showing 



the reliability of the method beyond the hopping expansion of eq. ITO 



At K = 0.16 we only show the full results for n?, = 1,2,3. The case 
Uf) = 2 shows a different shape with respect to the other two and indicates 
the relevance of mj2 terms. The simulations with Ub = 5 are only shown for 
the lowest values of p. We have observed that in this case the correlation 
times increase by almost two order of magnitude showing a severe critical 
slowing down. This may be ascribed to the delicate cancellation of "mag- 
netic" effects among three out the five bermion flavours. We plan to further 
study algorithms to improve the thermalization in a multiboson environment. 
In the region where we can perform the separate extrapolations of = 1,3 
and nb = 2, 5 results we observe again a reasonable consistency among them. 

We have shown that, by suitably modifying the boundary conditions, 
the bermion method can be applied to the study of full QCD at finite tem- 
perature and compared successfully with existing results. A detailed analysis 
of the nature of the phase transition, of the form of the effective action and 
of the possible interruption of the phase transition line for lower values of 
quark masses is currently under study [§]. 
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Table 1: The effective flavour numbers nf,mfi and m/2 obtained with bound- 
ary conditions Zj in theories with rib bermion fields. 
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Figure 1: The critical value of j3 in bermion simulations on a 16'^ x 2 lattice 
as a function of SrufiK^ is compared with the hopping expansion prediction 
[2] (see eq. 11). 
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Figure 2: The average space-space plaqucttc 1 — P^- and the absohite vahie 
of the Polyakov loop |L| on a 16^ x 4 lattice for k = 0.12 for four simulations 
with negative nj (corresponding, from right to left, to cases e, c, 6, a of table 
1) and the extrapolations to the fermion case n/ = 2 discussed in the text 
(★: extrapolation from cases a and c; x : extrapolation from cases h and e) . 
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Figure 3: The same as in fig. 2 for k = 0.14. 
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Figure 4: The same as in fig. 2 for k = 0.16. 
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Figure 5: The difference between space-space and space-time plaquettes for 
K = 0.14 for four simulations witli negative Uf (corresponding, from riglit 
to left, to cases e. c. b, a of table 1) and the extrapolations to the fermion 
case Hf = 2 discussed in the text extrapolation from cases a and c; x: 
extrapolation from cases b and e). 
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